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First evidence of a cyclotron feature in an anomalous X-ray pulsar 
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We report on the results of the longest uninterrupted observation of the Anomalous X-ray Pulsar 
IRXS J170849-400910 obtained with the BeppoSAX satellite in August 2001. The best fit phase-averaged spec- 
trum was found to be an absorbed power law plus blackbody model, with photon index P ~ 2.4 and a black body 
temperature of kTtt ^ 0.4 keV. We confirm the presence of significant spectral variations with the rotational phase 
of the pulsar. In the spectrum corresponding to the rising part of the pulse we found an absorption-like feature at 
~ 8.1 keV (significance level of 4a), most likely due to cyclotron resonant scattering. The centroid energy converts 
to a magnetic field of 9 x lO" G and 1.6 x lO^'^ G in the case of electrons and protons, respectively. If confirmed, 
this would be the first detection of a cyclotron feature in the spectrum of an anomalous X-ray pulsar. 



1. Introduction 

Anomalous X-ray pulsars (AXPs) are charac- 
terized by spin periods in the range of 5-12 s, 
steady spin down (~ 10~^^ ss~^), steep and soft 
X-ray spectra with luminosities exceding by sev- 
eral orders of magnitude their spin-down lumi- 
nosities (Mereghetti & Stella 1995). All five con- 
firmed AXPs lie in the "galactic" plane and two 
(or three), are associated with supernova rem- 
nants. AXPs show no evidence for a companion 
and are thus believed to be isolated neutron stars 
either having extremely strong magnetic dipole 
fields (~ 10^"* - IQi^G; "magnetars"; Duncan & 
Thompson 1992) or accreting from a residual disk 
(Alpar 2001). For recent rewiews see Mereghetti et 
al. 2002, and references therein. 

IRXS J170849-400910 was discovered with 
ROSAT (Voges et al. 1996) and ~ lis pulsations 
were found in its X-ray flux with ASCA (Sug- 
izaki et al. 1997). This sourse presents a strange 
timing behaviour, it is a quite stable rotator ex- 
cept few glitch type events that occurred in its 
secular spin-down (Kaspi et al. 2000, Dall'Osso 
et al. 2003, Kaspi & Gavriil 2003). Searches for 
an optical counterpart ruled out the presence of 
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a massive companion (Israel et al. 1999); an IR 
counterpart has been recently proposed (Israel et 
al. 2003). There is no evidence of pulsed radio 
emission from the AXP with an upper limit of 
70/iJy on the pulsation amplitude (Israel et al. 
2002). 

Here we report on the longest observation 
(200.0 ks) of this source, made with BeppoSAX 
in August 2001. Analysing spectral data we dis- 
covered the first absorption-like feature ever seen 
in an AXP, probably due to cyclotron resonant 
scattering. For more details about the observa- 
tion and the analysis see Rea et al. 2003. 

2. Timing and spectral results 

The timing analysis allowed us to estimate a 
spin period of 11.000563 ±0.000005 s. The source 
shows an energy-dependent profile (Fig. 3), in 
particular the pulse minimum shifts from a phase 
of ~0.0 in the lowest energy light curve (0.1- 
2keV) to ~0.3 in the 6-10 keV hght curve. Cor- 
respondingly, the pulsed fraction decreases from 
-30% to ~17%. 

The energy spectra were well fit with an ab- 
sorbed blackbody plus a power law model. The 
best fit of the phase-average spectrum gave a re- 
duced of 0-95 for 298 degree of freedom (dof) 
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Figure 1. Spectral parameters variability with 
phase, using an absorbed black body plus power- 
law model. Filled squares represent the spectral 
parameters after the addition a of cyclotron line 
model in the 0.4 - 0.58 phase interval. Note that 
the relevant points are slightly shifted in phase 
for clarity. 



for the following parameters: column density of 
Nh = (1.36±0.06) X 10^2 cm-2, a blackbody tem- 
perature of kTbb = 0.44 ± 0.01 kcV (blackbody ra- 
dius of Rbb = 6.6 ± 0.4 km, assuming a distance 
of 5 kpc) and a photon index of F = 2.40 ± 0.06 
(all error bars in the text arc 90% confidence). 
The unabsorbed flux in the 0.5-10 keV range was 
1.87x10"^° erg cm~2 s~^ corresponding to a lumi- 
nosity of 5.6 X 10'^^ erg s^^ (for a 5 kpc distance). 
In the 0.5-10 keV band the blackbody component 
accounts for ~ 30 % of the total unabsorbed flux 
(see Fig. 2). 



Figure 2. The first two panels show the pulsed 
fractions of the first and the second harmonics 
of the spin period in function of energy. In the 
third panel is presented the relative intensity of 
the fluxes of the power-law and the black body 
components and in the last one, of the power-law 
and the total flux, all in function of energy. 



In order to search for spectral features we made 
a pulse phase spectroscopy and accumulate spec- 
tra in six different phase intervals (see Fig. 1 and 
Tab. 1). A significant variation of the spectral 
parameters with pulse phase was clearly detect 
(especially for F) and confirm what previously 
reported in literature (Israel et al. 2001). In 
all intervals but one, an acceptable fit was ob- 
tained with the absorbed power law plus black 
body model (reduced in the 0.9-1.1 range); a 
reduced of 1.2 was instead obtained in the 0.4- 
0.58 phase interval. In the latter case, the data 
were systematically below the best fit model in 
the ~ 7.8-8.4 keV range (see Fig. 4b). 
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Figure 3. MECS light curves of 

IRXS J170849-400910 folded at the best 
spin period (two pulse cycles are shown) for six 
energy bands: (a) 0.1-2 keV; (b) 2 3 keV; (c) 
3-4 keV; (d) 4-5 keV; (e) 5-6 keV; (f) 6 10 keV. 
The vertical lines mark the phase interval in 
which the absorption-like feature was detected. 



We tried to fit three different models: a Gaus- 
sian, an absorption edge and a cyclotron feature. 



While the inclusion of a Gaussian or an absorp- 
tion edge did not lead to a significant improve- 
ment of the fit, the cyclotron model (CYCLABS 
in the XSPEC package, see Mihara et al. 1990 for 
details) led to a reduced of 1.09, corresponding 
to an F-tcst probability of 1.8 x 10~^ implying to 
a single trial significance of 4.5cr or 4a after cor- 
rection for the six spectra that we analysed (see 
Fig. 4 and Tab 1). Moreover, taking into account 
also the possibility of find the feature in all 1-10 
keV range, we found a 3.5cr confidence level. 

3. Discussion 

During the analysis of the longest BeppoSAK 
observation of IRXS J170849-400910 we discov- 
ered an absorption-like feature at an energy of 
^8.1 keV in a pulse phase interval corresponding 
to the rising part of the ~ 11 s pulse. This feature 
was best fit by a resonant cyclotron feature model 
with a centroid energy of ~8.1keV and an equiv- 
alent width of ~460eV (see Fig. 4 and Tab. 1). 

The detection of an RCF in a specific pulse- 
phase interval and superposed to an X-ray contin- 
uum that varies with the pulse phase is reminis- 
cent of the behaviour seen in standard accreting 
pulsars in X-ray binaries (Wheaton et al. 1979). 
If interpreted as an electron resonant feature at 
the base of the accretion column, the feature at 
~ 8.1 keV implies a neutron star surface magnetic 
field of ~ 9.2 X 10^^ Gauss (using a gravitational 
redshift z=0.3) . This value is just slightly lower 
than that measured for electron RCFs in typical 
accreting X-ray pulsars (see Fig. 5); more inter- 
estingly it is close to that required by models for 
AXPs which involve residual disk accretion in the 
spin-down regime. In this context, one can solve 
the torque equation (see e.g. Eq. 11.35 in Hen- 
richs 1983) by exploiting the measured value of 
P and range of accretion luminosity derived from 
plausible distances (5-10 kpc). The surface mag- 
netic (dipolc) fi(4(l ol)tained in this way is 0.6- 
1.1x10^^ G (corrisponding to a fastness param- 
eter range of Ug = 0.57 — 0.54, a typical value 
for the spin-down accretion regime; see Ghosh & 
Lamb 1979 and Henrichs 1983). 

The agreement of this estimate with the mag- 
netic field inferred from the electron RCF inter- 
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pretation is intriguing, especially in consideration 
of the other analogies with the pulse-phase spec- 
tral dependence of conventional accreting X-ray 
pulsars. By contrast, if an electron RCF arose 
somehow at the polar caps of a rotation powered 
pulsar, a B field strength of 9.2xlO^^G would 
be in the range of many radio pulsars and yet 
much lower than that required to spin-down at 
the observed rate through magnetic dipole radi- 
ation (~ 5 X 10^^ G, indeed this was one of the 
motivations for magnetar model, see below). 

There is a clear correlation between the width 
and centroid energy of the electron RCFs in ac- 
creting X-ray pulsars (see Fig. 5, extending the 
results of Orlandini & Dal Fiumc 2001). The val- 
ues from IRXS J170849-400910 and SGR1806- 
20 are in good agreement with such a relation. 
The modest range of width to centroid energy ra- 
tio implied by this indicates that magnetic field 
geometry effects at the neutron star surface likely 
dominate the RCF width (on the contrary tem- 
perature and particle mass would alter this ra- 
tio). This, in turn, suggests that similar (rel- 
ative) ranges of surface magnetic field strength 
are "sampled" by RCFs in accreting X-ray pul- 
sar and, by extension, RCFs in AXPs and SGRs. 

Alternatively the RCF might be due to pro- 
tons. For the magnetic field strengths forseen 
in the "magnetar" scenario, proton cyclotron fea- 
tures (if any) are expected to lie in the classical 
X-ray band (0.1-10 keV; Zane et al. 2001; Lai & 
Ho 2002). A proton RCF feature at ~ 8.1 keV 
would correspond to surface field of 1.6xlO"'^^G 
(z=0.3). The fact that this value is ~ 3 times 
higher than the surface field derived from the 
usual magnetic dipole spin -down formula should 
not be of concern. According to the magnetar 
model, the magnetic field at the star surface and 
its vicinity is dominated by higher order multipole 
field components. At large radii the dipole com- 
ponent, responsible for the secular spin-down, 
dominates. It is thus expected that a proton 
RCF feature, sampling the (total) surface mag- 
netic field strength, provides a higher value than 
the mere dipole component. 

Other interpretations of the feature at ~ 
8.1 keV appear less likely. Firstly, fitting an edge 
or line due to photo-electric absorption provides 
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Figure 4. MFCS and LEGS spectra from the 0.4- 
0.58 phase interval fit with the "standard model" 
(the sum of a blackbody and power law with ab- 
sorption) plus a cyclotron line. Residuals are rela- 
tive to the "standard model' alone in order to em- 
phasize the absorption-like feature at ~ 8.1 keV: 

(a) the BeppoSAX observations merged together; 

(b) the 2001 observation alone; (c) the phase in- 
tervals contiguous to that showing the cyclotron 
absorption feature in the merged observations. 



a less pronounced improvement of the fit than 
the RCF model. Secondly, an edge by iron at a 
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Table 1 



Spectral parameter values of the phase resolved spectroscopy 



Phase 


0.0-0.26 


0.26-0.4 


0.4-0.58 


0.58-0.7 


0.7-0.84 


0.84- -1.0 


iVH(xlO^^ cm"^) 


1.26 ±0.11 


1.6±0.1 


1.38 ±0.12 


1.24 ±0.14 


1.35 ±0.12 


1.1±0.1 


F 


2.2 ±0.1 


2.9 ±0.1 


2.62 ±0.13 


2.60 ±0.13 


2.30 ±0.15 


2.03 ±0.06 


kTbb (keV) 


0.433 ±0.017 


0.424 ±0.013 


0.46 ±0.01 


0.49 ± 0.03 


0.465 ±0.018 


0.46 ±0.02 


Rbb (d = 5 kpc; km) 


6.0 ±0.3 


7.2 ±0.4 


6.7 ±0.4 


6.1 ±0.3 


6.9 ±0.3 


5.0 ±0.3 


Ecyclabs (keV) 






8.1 ±0.1 








Line Width (keV) 






0.2 ±0.1 








Line Depth ( keV) 






0.8 ±0.4 








xVd.o.f. 


1.08 


0.98 


1.09 


0.99 


1.09 


1.05 



sufficiently large distance from the neutron star 
that energy shifts are negligible would require a 
high overabundance of this element and interme- 
diate ionisation stages (such as C-like iron). Yet 
it has long been known that the photoionisation 
equilibrium of such a plasma is unstable (Kro- 
lik & Kallman 1984; Nagase 1989). The energy 
of an ion feature forming in the neutron star at- 
mosphere would be drastically altered by mag- 
netic field effects (see Mori & Hailey 2002 and 
references therein). In this and the above inter- 
pretations, however, it would also be difficult to 
explain why an ion feature is observed only over 
a restricted range of pulse phases. 

4. Conclusion 

In conclusion, we found an absorption-like 
feature in the BeppoSAX X-ray spectrum of 
IRXS J170849-400910 taken during the rising 
phase of the ~lls pulse, which is best fit by a 
cyclotron resonant scattering model. If this inter- 
pretation is correct, the centroid energy translates 
into a model-indipendent magnetic field strength 
of ~1.6xl0i^G or ~9.2xlO"G depending on 
whether protons or electrons, respectively, are re- 
sponsible for the feature. 

Further observations and polarization measure- 
ments might choose beetwen the two case. If 
the line is due to protons this result will give 
the confirmation that AXPs are magnetars, if in- 
stead, is due to electrons, this depress the mag- 
netar model in favor of accretion-driven models 
for analogy with other accreeting pulsar magnetic 
field strenghts known so far. 
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Figure 5. Line width vs. centroid en- 
ergy from a sample of accreting X-ray pulsars 
with electron RCFs (Orlandini & Dal Fiume 
2001), 1RXSJ170849-400910 (this Letter) and 
SGR 1806-20 (Ibrahim et al. 2002). 



